R代码:GEO单个样本文件的合并 |
您所在的位置:网站首页 › merge data file › R代码:GEO单个样本文件的合并 |
在GEO中下载MINiML formatted family file(s) setwd("C:/AW") PROBE=c("26540_at","76571_s_at")#目标基因的两个探针 PROBE=as.data.frame(PROBE) colnames(PROBE)="probe" a=list.files("GSE39582_family") #提取GSE39582_family文件夹下所有文件名 n=length(a) #GSE39582_family文件夹文件数量 for(i in 1:n){ b[i]=substring(a[i],1,10) #在文件名中取前10即样本名 } dir=paste("./GSE39582_family",a,sep="/") merge.data=read.delim(file=dir[1],header=F) #提取第一个文件 colnames(merge.data)=c("probe",b[1]) #列名 for(i in 2:n){ new.data=read.delim(file=dir[i],header=F) colnames(new.data)=c("probe",b[i]) merge.data=merge(merge.data,new.data,by="probe") #根据“probe”列合并文件 } write.csv(merge.data,file = "./GSE39582_family/merge.csv",row.names=F) #输出合并文件 expression=merge(merge.data,PROBE,by="probe",all=F) #提取合并文件中目标基因列 write.csv(expression,file = "./GSE39582_family/GSE39582C1.csv",row.names=F) #输出目标基因文件 参考:https://blog.csdn.net/weixin_40645816/article/details/82110057 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |